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In this paper we study the sensitivity of the Large Sky Area Multi-Object Fiber Spectroscopic 
Telescope (LAMOST) project to the determination of cosmological parameters, employing the Monte 
Carlo Markov Chains (MCMC) method. For comparison, we first analyze the constraints on cosmo- 
logical parameters from current observational data, including WMAP, SDSS and SN la. We then 
simulate the 3D matter power spectrum data expected from LAMOST, together with the simu- 
00 "■ lated CMB data for PLANCK and the SN la from 5-year Supernovae Legacy Survey (SNLS). With 

the simulated data, we investigate the future improvement on cosmological parameter constraints, 
emphasizing the role of LAMOST. Our results show the potential of LAMOST in probing for the 
£NJ ■ cosmological parameters, especially in constraining the equation-of-state (EoS) of the dark energy 

and the neutrino mass. 
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INTRODUCTION 



The measurement of the large scale galaxy clustering has been an important probe in constraining the cosmological 
models. The large scale structure (LSS) measurements have made remarkable progresses by the observational efforts 
such as 2dFGRS and Sloan Digital Sky Survey (SDSS), which have provided an accurate measurement of the galaxy 
I ■ power spectrum and given a robust constraint on cosmological parameters [l], Q • 

The LAMOST [| project is a 4m quasi-meridian reflecting Schmidt telescope laid down on the ground. It has a 5 
■ degree field of view, and may accommodate as many as 4000 optical fibers and the light from 4000 celestial objects will 
be led into a number of spectrographs simultaneously. Thus the telescope will be the one that possesses the highest 
l—— spectrum acquiring rate in the world. The spectroscopic survey which contains the information about the radial 
| positions of galaxies, can probe the 3D distribution of galaxies effectively. In this paper, we study the sensitivity 
of LAMOST to the determination of cosmological parameters with the simulated galaxy power spectrum. In our 
£S) ' analysis, we also consider the simulated observations for the future CMB and SN la measurements from PLANCK 
and the 5-year SNLS, which are presumably conducted during the same time period as LAMOST survey. Our results 
indicate that the LAMOST has the promising potential in probing for the cosmological parameters, especially in 
. ■ constraining on the EoS of the dark energy and the neutrino mass. 

The paper is organized as follows: In Section II we describe the method of fitting and the simulation technique. In 
. section III, we present the results and discussions. The last section contains a summary. 
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METHODOLOGY 



In this section, we introduce the method and the fitting procedure. For the dynamical dark energy model, we 



choose the parametrization given by [4|: 

w(a) = wo + w a (l - a) , (1) 

where a — 1/(1 + z) is the scale factor and w a = —dw/da characterizes the "running" of the EoS (Run w henceforth). 
For the ACDM model, wq = — 1 and w a = 0. 

When using the MCMC global fitting strategy to constrain cosmological parameters, dark energy perturbations 
should be taken into account properly, especially for models with time evolving EoS of dark energy. This issue has 
been realized by many researchers including the WMAP group @, 0, H|. However, when the parameterized EoS 
crosses —1, one cannot handle the dark energy perturbations based on quintessence, phantom, k-essence and other 
non-crossing models. By virtue of quintom the perturbations at the crossing points are continuous. Thus we have 
proposed a technique to treat dark energy perturbations in the whole parameter space. 

In this study, we have modified the publicly available Markov Chain Monte Carlo package CosmoMC[10( to include 
the dark energy perturbations. For handling the parametrization of the EOS getting across -1, firstly we introduce a 
small positive constant e to divide the full range of the allowed value of the EOS w into three parts: 1) w > — 1 + e; 
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2) — 1 + e > w > — 1 — e; and 3) w < — 1 — e. Working in the conformal Newtonian gauge, the perturbations of DE 
can be described by 

5 = -(l+w)(9-3$)-3H(c 2 s -w)5 , (2) 

7/'; Ft 

9 = -H(l-3w)6-- 9 + k 2 (-^— + *) . (3) 

v ' l + w A + w ' w 

Neglecting the entropy perturbation, for the regions 1) and 3), the EOS does not across —1 and the perturbation 
is well defined by solving Eqs. (|2l3|) . For the case 2), the perturbation of energy density 5 and divergence of velocity, 
9, and the derivatives of 5 and 9 are finite and continuous for the realistic quintom DE models. However for the 
perturbations of the parameterizations, there is clearly a divergence. In our study for such a regime, we match the 
perturbations in region 2) to the regions 1) and 3) at the boundary and set 

(5 = , 9 = 0. (4) 

In our numerical calculations we limit the range to be \Aw = e| < and find our method to be a very good 

approximation to the multi- field quintom. More detailed treatments can be found in Ref. 0, |6[ . 

Furthermore, we assume purely adiabatic initial conditions and a flat universe. The parameter space we begin with 
for the numerical calculation is: 

P = (cj b ,w c ,e s ,T, u> , w a ,n s ,hi(lO w A s )) , (5) 

where u>b = Vt^h 2 and u c = fl c h 2 with fib and fl c being the baryon and cold dark matter densities relative to the 
critical density, respectively, <d s is the ratio (multiplied by 100) of the sound horizon at decoupling to the angular 
diameter distance to the last scattering surface, and r is the optical depth. In Eq.([5]), A s and n s characterize the power 
spectrum of primordial scalar perturbations. For the pivot scale of the primordial spectrum we set fc* = 0.05Mpc _1 . 

In our calculations, we take the total likelihood to be the products of the separate likelihoods (d) of CMB, LSS 
and SNIa. Defining \\ « = — 21og£.;, we then have 

2 _ 2 i 2 i 2 /f.\ 

AL,total — AL,CMB t XL,LSS t XL,SNIa ■ \°) 

If the likelihood function is Gaussian, \\ coincides with the usual definition of x 2 U P to an additive constant corre- 
sponding to the logarithm of the normalization factor of C. 

The data used for current constraints include the three-year WMAP (WMAP3) 1 Temperature- Temperature (TT) 
and Temperature-Polarization (TE) power spectrum 0, HE! EH as well as the smaller scale experiments, including 



Boomerang-2K2 [T3], CBI [3], VSA and ACBAR [2c|, the SDSS luminous red galaxy (LRG) sample Q and 
2dFGRS 0], and recently released ESSENCE (192 sample) data 22|, [23[. For the LSS power spectrum, we only use 
the data in the linear regime up to k ~ QAhMpc^ 1 . In the calculation of the likelihood from SNIa we have marginalized 
over the nuisance parameter [2l|. Furthermore, we make use of the Hubble Space Telescope (HST) measurement of 
the Hubble parameter Hq = lOOh km s _1 Mpc -1 by multiplying the likelihood by a Gaussian likelihood function 
centered around h = 0.72 and with a standard deviation a = 0.08. We also impose a weak Gaussian prior on the 



baryon density f^ft 2 = 0.022 ± 0.002 (1 a) from the Big Bang Nucleosynthesis [25|, and a cosmic age tophat prior as 
10 Gyr < t < 20 Gyr. 

For the future data, we consider the measurements of LSS from LAMOST, the CMB from PLANCK (HI and the 



SN la from 5-year SNLS[27| 



For the simulation of LAMOST, we mainly simulate the galaxy power spectrum. We consider two sources of statical 
errors on the power spectrum measurements: the sample variance and the shot noises which due to the limited number 
of independent wavenumber sampled from a finite survey volume and the imperfect sampling of fluctuations by the 
finite number of the galaxies respectively, [28| 

where V is the survey volume and n is the mean galaxy density. From the more conservative estimation, we know 
that the redshift distribution of main sample of LAMOST is between and 0.6 and the mean redshift is around 0.2. 



1 With the newly released 5-year WMAP datallll. ll2H . we have checked that the new data will not significantly change the results. 
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TABLE I. Constraints on cosmological parameters from the current observations and the future simulations. For the current constraints 
we have shown the mean values la (Mean) and the best fit results together. For the future mocked data we list the standard deviation 
(SD) of these parameters with fiducial model I (FMI) and fiducial model II (FMII). In order to highlight the contribution from 

LAMOST, we compare the results with/without LAMOST. 





Current for ACDM 


Future (SD with FMI) 


Current for Run w 


Future (SD with FMII) 




Best Fit 


Mean 


PLANCK+ SNLS 


PLANCK+ SNLS 
+LAMOST 


Best Fit 


Mean 


PLANCK+SNLS 


PLANCK+ SNLS 
+LAMOST 


wo 


-1 


-1 


0.118 


0.100 


-1.16 


-1 03+° 15 


0.0899 


0.0598 


W a 








0.522 


0.417 


0.968 




0.332 


0.162 


n A 


0.760 


762+ - 015 


0.0115 


0.00547 


0.756 


760+ - 017 
U -' DU -0.018 


0.0125 


0.00460 




73.1 


73.31 \- 6 r 


1.594 


0.828 


70.3 


71.2j£3 


1.840 


0.673 




0.769 


0.755 ± 0.031 


0.0223 


0.0174 


0.634 


0.675 ± 0.068 


0.0299 


0.0220 



So for simplicity, in our study, we simulate the power spectrum of the galaxies at z = 0.2 that can be got from these 
galaxies. The survey area is 15000 deg 2 and the total number of galaxies within the survey volume is 10 7 [3j. We only 
consider the linear regime, namely the maximum k we consider is k ~ 0.1 h Mpc" 1 . As we know that, the galaxy 
power spectrum P(k) in EQ. ([7]) is 

P(k) = b 2 Pm (k), (8) 

where p m {k) is the linear matter power spectrum, and here we take b as a constant 6=1 when simulating the data 
and when using the galaxy power spectrum to constrain cosmological parameters, we take b as a free parameter and 
marginalize over it. 

For the simulation with PLANCK, we follow the method given in our previous paper (29[. We mock the CMB TT, 
EE and TE power spectrum by assuming the certain fiducial cosmological model. For the detailed techniques, please 
see our previous companion paper (29[. We have also simulated 500 SN la according to the forecast distribution of 
the SNLS [30]. For the error, we follow the Ref.j3l| which takes the magnitude dispersion 0.15 and the systematic 
0.02 x z/1.7. The whole error for each data is given by: 



/ 0.15 2 

Vmaga(Zi) = \ V 2 ys( z i) + ~ > (9) 

y Tli 

where rii is the number of supernova of the i'th redshift bin. 

As pointed out in our previous works H, 33, 34 1, the cosmological parameters are highly affected by the dark energy 



models due to the degeneracies among the EoS of DE and other parameters. Therefore, in our study of this paper, 
we choose two fiducial models with different dark energy properties: ACDM model (fiducial model I henceforth) and 
dynamical dark energy (fiducial model II henceforth) with time evolving EoS. The parameters of the two sets of 
fiducial models are obtained from the current observational data. 



RESULTS AND DISCUSSIONS 



In Table I, we present the numerical results of the constraints on the cosmological parameters from the current 
data and the error forecast from the simulated data. To show the importance of LAMOST, we compare the two sets 
of results, one from PLANCK + SNLS, the other from PLANCK + SNLS + LAMOST. As we know, the matter 
power spectrum is directly related to the horizon size at matter-radiation equality, in turn the matter power spectrum 
will make accurate measurement of Sl m h. On the other hand, there are degeneracies between £l m h and the other 
cosmological parameters, e.g. Qa, Hq, wq, w a and so on, hence the tight constraint on Q m h will be helpful for breaking 
these degeneracies and improve the constraints on these cosmological parameters. For example, in Table I, one can 
find the constraints on and H are tightened a lot by including LAMOST. Also from Figure 1, one can see the 
constraits on the age of universe is also shrank obviously, this is because the age is directly related to the hubble 
constant and f2 m . 

In figure [TJ we plot the 2— D cross correlation and 1— D probability distribution of some of the basic cosmological 
parameters in Eq.([5]) and also some of the reduced parameters. The black solid lines are given by fitting with the 
simulated PLANCK and SNLS, and the red solid lines are provided by including the simulated LAMOST data. From 
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TABLE II. Constraints on neutrino mass from the current observations and the future simulations. We have shown the 2<r 
upper limits. In order to highlight the contribution from LAMOST, we compare the results with/without LAMOST. 



Current 


Future 








PLANCK+SNLS 


PLANCK+SNLS+LAMOST 


ACDM 


<0.958eV 


FMI 


<0.957eV 


<0.377eV 


Run w 


<1.59eV 


FMII 


<0.915eV 


<0.346eV 



the comparison, we find that LAMOST have promising potentials in constraining cosmological parameters, such as, 
the EoS of dark energy, the dark energy density budget the age of universe, erg and the Hubble constant. 

In order to see explicitly the effect of LAMOST on dark energy constraints, in figure [2j we plot the 2a confidence 
level contours on wq and w a . The black solid line is given by the current constraints and the red solid line is given 
by fitting with the simulated data of PLANCK and SNLS 5 year data with fiducial model II, while the red dashed 
line is given by including the simulated LAMOST data. This comparison shows clearly LAMOST will contribute 
significantly in tightening the constraints on the EoS of dark energy. Numerically we find the best fit model with 
the current data is given by the dynamical quintom model with EoS across -1, however the cosmological constant is 
within I a confidence level. The future PLANCK measurement and 5-year SNLS SN la will be able to distinguish the 
cosmological constant from the dynamical model at 2a confidence level, while LAMOST can improve this sensitivity 
significantly at 3.3<r. The blue solid lines and blue dashed lines show the comparison between the results with and 
without LAMOST for fiducial model I. 

On the other hand, for the parameter related to the inflation models, such as n s and A s , the constraints are mainly 
from PLANCK, as pointed out in our previous paper [2{|. Adding in LAMOST can not further tighten the constraints. 
We have also done another analysis with the additional parameters a s and r, and obtained the similar conclusion, 
where a s characterizes the running of the primordial power spectrum index and r is the ratio of tensor to scalar 
perturbations. 

Now we study the cosmological constraint on the neutrino mass by adding in a new parameter /„ in Eq.([5]). The 
parameter /„ is the dark matter neutrino fraction at present, namely, 

f =-£^ = ^ do) 

Pdm 93.105 eV fl c h 2 ' v ; 

where £to„ is the sum of the neutrino masses. In this study, the mocked data we use are generated by assuming the 
massless neutrino i.e. /„ = in the fiducial models. Consequently the constraints on /„ should be regarded as the 
upper limits of the neutrino mass which the future observations will be sensitive to. It is well known that the massive 
neutrinos modify the shape and amplitude of the matter power spectrum, and also the epoch of matter-radiation 
equality, angular diameter distance to the last scattering surface. Thus they leave imprints on the observations of 
CMB and LSS. In Table II, we provide the constraints on neutrino mass from the current observations and the future 
simulated data. For the current data 2 , within the framework of the ACDM model, we get Em„ < 0.958 eV (95%) 
which is consistent with the result in Ref. [l[ . For the time evolving EoS of dark energy model, this limit is relaxed to 
Y>m v < 1.59 eV (95%), due to the degeneracy between the dark energy parameters and the neutrino mass(33l. l36j|. 

With the simulated data, in figure [3l we illustrate the one dimensional probability distribution of the total neutrino 
mass ^2m u . The black solid line is given by the current constraints, the red solid line is given by fitting with the 
simulated PLANCK + SNLS with fiducial model II and the red dashed line is given by including the simulated 
LAMOST. The blue solid line and blue dashed line are the results obtained with the fiducial model I. Our results 
show that the LAMOST can provide a more stringent constraint on the neutrino mass. For example, the 2a neutrino 
mass limit is changed from 0.957eV to 0.377ey by including the simulated LAMOST data with the fiducial model I. 



2 Usually Lyman-a data will give stringent constraint on the neutrino mass, however, its systematics are quite unclear currently. In our 
global analysis to be conservative, we have not included it [12ll35l|. 
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FIG. 1: one-dimensional distributions and two-dimensional 68% and 95% limits on the cosmological parameters. The black 
solid lines are obtained with the simulated PLANCK + SN la and the red solid lines are from PLANCK + SN la + LAMOST. 



SUMMARY 



In this paper we have studied the sensitivity of LAMOST project to the determination of the cosmological param- 
eters. With the simulated 3D matter power spectrum of LAMOST, in combination with the future PLANCK data 
and 5-year SNLS data, we have obtained the constraints on the various parameters by employing the MCMC method. 
Our results show the potential for LAMOST in constraining the cosmological parameters, especially on the EoS of 
dark energy and the neutrino mass. 

We have performed our analysis in flat universe, however, if we take the curvature ilk into consideration, namely if 
we free ilk in the global analysis, the basic conclusion will not change. That is to say, we can also find the potential of 
the future LAMOST data in determining cosmological parameters, however, the specific contours of each cosmological 
parameters will be enlarged for the additional degree of freedom, and more relevant discussion can be seen in our 
previous paper 34], in which we have implemented the global fitting with the observational data for non-flat universe. 
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FIG. 2: 2-d joint 68% and 95% confidence regions of the parameters tun and w a for flat universe. The black solid line is given 
by the current constraints, the red solid line comes from the simulated data of PLANCK and 5year SNLS with fiducial model 
II, while the red dashed line is by combining the simulated LAMOST data. The blue solid line and blue dashed line are the 
results with the fiducial model I with/without LAMOST respectively. 




FIG. 3: one-dimensional constraints on the neutrino mass. The black solid line is given by the current data, the red solid line 
is given by fitting with the simulated PLANCK + SNLS with fiducial model II and the red dashed line is given by including 
the simulated LAMOST. The blue solid line and blue dashed line are the results obtained with the fiducial model I. 
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